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Abstract 

We report new radial velocities of the TrES-4 transiting planetary system, including observations of 
a full transit, with the High Dispersion Spectrograph of the Subaru 8.2ni telescope. Modeling of the 
Rossitcr-McLaughlin effect indicates that TrES-4b has closely aligned orbital and stellar spin axes, with 
A = 6.3° ±4.7°. The close spin-orbit alignment angle of TrES-4b seems to argue against a migration history 
involving planet-planet scattering or Kozai cycles, although there are two nearby faint stars that could be 
binary companion candidates. Comparison of our out-of-transit data from 4 different runs suggest that 
the star exhibits radial velocity variability of ~20 m s~^ in excess of a single Keplerian orbit. Although 
the cause of the excess radial velocity variability is unknown, we discuss various possibilities including 
systematic measurement errors, starspots or other intrinsic motions, and additional companions besides 
the transiting planet. 

Key words: stars: planetary systems: individual (TrES-4) — stars: rotation — stars: binaries: 
general — techniques: radial velocities — techniques: spectroscopic 



1. Introduction 

Transiting planets provide us with valuable opportu- 
nities to characterize exoplanetary systems. One such 
opportunity is to measure the Rossiter-McLaughlin ef- 
fect (hereafter the RM effect: Rossiter 1924, McLaughlin 
1924), an apparent radial velocity anomaly during a plan- 
etary transit, which is caused by a partial eclipse of the 
rotating surface of the host star. By measuring and mod- 
eling this effect, one can measure the sky-projected an- 
gle between the stellar spin axis and the planetary or- 
bital axis. Many theoretical investigations of the RM ef- 
fect have been presented (e.g., Ohta et al. 2005; Gimenez 
2006; Gaudi & Winn 2007; Hirano et al. 2010), and obser- 
vations of the RM effect have been reported for about 20 
transiting planetary systems (for the most recent compila- 
tion, see Jenkins et al. 2010 and references therein). One 
of the main theoretical motivations to observe the RM ef- 
fect is that the observed degree of spin-orbit alignment is 
thought to be connected with the migration history of the 
transiting planet. 

The most frequently discussed planetary migration 
mechanisms are (1) gravitational interaction between a 
protoplanetary disk and a growing planet (disk-planet in- 
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teraction models, e.g., Lin & Papaloizou 1985; Lin et al. 
1996; Ida & Lin 2004), (2) gravitational planet-planet 
scattering and subsequent tidal evolution (planet-planet 
scattering models, e.g., Rasio & Ford 1996; Marzari & 
Weidenschilling 2002; Nagasawa et al. 2008; Chatterjee 
et al. 2008), or (3) the Kozai mechanism caused by a dis- 
tant massive companion and subsequent tidal evolution 
(Kozai migration models, e.g., Wu & Murray 2003; Takeda 
& Rasio 2005; Fabrycky & Tremaine 2007; Wu et al. 2007). 
These scenarios are not necessarily mutually exclusive, 
but to the extent that they are, disk-planet interaction 
models would predict small orbital eccentricities and good 
spin-orbit alignments, while planet-planet scattering mod- 
els and Kozai migration models predict a broader range 
of eccentricities and spin-orbit alignment angles. Until 
about a year ago, all of the measurements indicated close 
alignments, but recently 6 transiting planets have been 
reported with significant misalignments: X0-3b (Hebrard 
et al. 2008; Winn et al. 2009b), HD 80606b (Moutou et al. 
2009; Pont et al. 2009a; Winn et al. 2009c), WASP-14b 
(Johnson et al. 2009; Joshi et al. 2009), HAT-P-7b (Narita 
et al. 2009a; Winn et al. 2009a), CoRoT-lb (Pont et al. 
2009b), and WASP-17b (Anderson et al. 2010). With this 
increase in the number of measurements, and the diversity 
of results, we are approaching the time when we may test 
the validity and applicability of the different planetary 
migration models. 
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The main target of this paper is TrES-4b, which is 
a transiting exoplanet discovered by Mandushev et ah 
(2007) (hereafter M07) in the course of the TrES survey, 
supplemented by Keck radial velocity measurements. The 
planet orbits an F8 (Daemgen et al. 2009) host star every 
3.55 days and is one of the most "inflated" hot Jupiters 
with a radius of about 1.8 i?jup, which places this planet 
to be one of the least density exoplanets ever discovered. 
Refined spectroscopic and photometric characteristics of 
the host star TrES-4 were presented by Sozzetti et al. 
(2009) (hereafter S09). The amplitude of the RM effect 
for TrES-4b was expected to be large, because of the large 
projected equatorial rotational velocity of the TrES-4 star 
{Vsmis = 8.5 km s~^; S09). In addition, Daemgen et al. 
(2009) have recently reported a possible companion star 
around the TrES-4 system. Although it has not yet been 
confirmed that the companion is a true physical binary as 
opposed to a chance alignment, a companion star would 
raise the possibility of migration via Kozai cycles, lending 
additional motivation to the study of the RM effect in this 
system. We note that we adopt "TrES-4" as the host star 
name and "TrES-4b" as the planet name in this paper, al- 
though the planet was originally named "TrES-4" by the 
discoverers (see MOT). The reason of our choice is because 
recent papers on this system (e.g., Daemgen et al. 2009, 
which we referred in this paper) often describe the host 
star as "TrES-4" . Thus we consider that it would be con- 
fusing for readers if we describe the planet as "TrES-4" in 
our paper, and we hope to avoid such confusions. 

In this paper, we present new measurements of the ra- 
dial velocity of TrES-4 made with the Subaru 8.2m tele- 
scope. Although TrES-4 is relatively faint {V = 11.6), 
the large aperture of the Subaru telescope has enabled us 
to measure radial velocities of TrES-4 with high precision. 
Our radial velocity dataset consists of 23 samples covering 
a full transit, and 8 samples obtained outside of transits on 
3 different nights. In addition to reporting the spin-orbit 
alignment angle of TrES-4b, we also report the observa- 
tion of radial velocity variation in excess of the previously 
observed single Keplerian orbit, and we confirm through 
direct imaging with the HDS slit viewer that there are 
candidate companion stars. 

The rest of this paper is organized as follows. Section 2 
summarizes our Subaru observations, and section 3 de- 
scribes analysis procedures of the RM effect of TrES-4b. 
Section 4 presents our main result on the spin-orbit align- 
ment angle of TrES-4b, and section 5 discusses possible 
causes of the additional radial velocity variation in this 
system. Finally, section 6 summarizes the findings of this 
paper. 

2. Observations 

We observed a full transit of TrES-4b with the High 
Dispersion Spectrograph (HDS: Noguchi et al. 2002) on 
the Subaru 8.2m telescope on UT 2007 July 13. In ad- 
dition, we measured radial velocities on UT 2007 August 
5, UT 2008 March 9, and UT 2008 May 30, when tran- 
sits were not occurring. For all observations, we employed 



the same setup. We used the standard I2a setup of the 
HDS, a slit width of 0"8 corresponding to a spectral res- 
olution of about 45,000, and the iodine gas absorption 
cell for precise differential radial velocity measurements. 
The exposure times for the radial velocity measurements 
were 12-30 minutes, yielding a typical signal-to-noise ratio 
(SNR) of approximately 80-100 per pixel. We processed 
the observed frames with standard IRAF^ procedures and 
extracted one-dimensional spectra. We note that in about 
half of the exposures, light from a nearby companion star 
(sec section 3.1) was also admitted through the slit. For 
this reason, when extracting the spectra we limited the 
aperture width of echelle orders to 3-4 pixels, to avoid 
including any significant fiux from the companion star. 
We computed relative radial velocities and uncertainties 
following the algorithm of Butler et al. (1996) and Sato 
et al. (2002), as described in Narita et al. (2007). We 
estimated the internal measurement uncertainty of each 
radial velocity based on the scatter of radial velocity so- 
lutions for ~4 A segments of each spectrum. The typical 
internal uncertainties were 15-20 m s~^, which is worse 
than typical cases having similar SNR due to relatively 
rapid rotation of TrES-4. The radial velocities and uncer- 
tainties are summarized in Table 1. 

3. Analyses 

3.1. HDS Slit Viewer Images 

The slit viewer of the Subaru HDS has a 512x512 CCD, 
providing unfiltered 60" x60" field of view images. During 
the HDS observations, we found two nearby companion 
stars in HDS slit viewer images. One is the companion star 
which Daemgen et al. (2009) also reported to the north of 
TrES-4, and the other is a newly-discovered star located 
in west-southwest. Figure 1 shows 4 magnified portions 
of slit viewer images obtained on UT 2007 July 13 (upper 
left), UT 2007 August 5 (lower left), UT 2008 March 9 
(upper right), and UT 2009 July 12 (lower right). North 
is up and east is left for these images, and the field of view 
is 20" X 20" . Our ability to perform astrometry and pho- 
tometry on these images is limited, since TrES-4 located 
on the slit. For the north companion star, we roughly esti- 
mate the separation angle as ^ l."56 (~ 13 pixels) and the 
position angle as ^ 0° (almost true north) . These findings 
are consistent with those of Daemgen et al. (2009). While 
the west-southwest companion star is located at the sepa- 
ration angle of ~ 7."8 (~ 65 pixels) and the position angle 
of ^ 249°. We also estimate that the companion stars 
are at least 4 magnitudes fainter than TrES-4 in visible 
wavelength. Although neither our data nor the data of 
Daemgen et al. (2009) can be used to tell whether or not 
the companion star is physically associated with TrES-4, 
it may be possible to do so with multiband IR (for exam- 
ple, JHK band) photometry with adaptive optics, which 

^ The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the U.S. National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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would be very useful to study common proper motion, 
color-magnitude relation, and spectral type of the com- 
panion stars. Such follow-up observations will enable us 
to investigate the binarity of the companion stars. 

3.2. Simulated Formula for the Rossiter-McLaughlin 
Effect 

We modeled the RM effect of TrES-4 following the 
procedure of Winn et al. (2005), as adapted for HDS 
by Narita et al. (2009a, 2009b) and discussed further 
by Hirano et al. (2010). We first created a synthetic 
template spectrum which matches the stellar properties 
of TrES-4 described by S09, using a model by Coelho 
et al. (2005). To simulate the disk-integrated spectrum 
of TrES-4, we applied a rotational broadening kernel of 
V&vo-Is = 8.5 km s~^ (S09) and adopted the quadratic 
limb-darkening parameters for the spectroscopic band as 
ui = 0.46 and U2 = 0.31 based on the tables of Claret 
(2004). To simulate in-transit spectra of TrES-4, we sub- 
tracted a scaled-down and velocity-shifted copy of the 
original unbroadened spectrum, which represents the hid- 
den part of the stellar surface. We created a collection 
of such simulated spectra for various values of the scal- 
ing factor / and the velocity-shift Wp, and computed the 
apparent radial velocity At; for each spectrum. We fitted 
Aw in (/, Wp) space and determined an empirical formula 
of the RM effect for TrES-4 as 



Aw — —fvp 



1.623-0.885 



y sin/. 



(1) 



3.3. Radial Velocity Modeling 



Since we did not have good transit light curves for TrES- 
4b. we fixed stellar and planetary parameters of TrES- 
4 to the values reported by S09 as follows; the stellar 
mass Me = 1.404 [Mq], the stellar radius Rs = 1.846 [i?©], 
the radius ratio Rp/Rg = 0.09921, the orbital inclination 
i = 82.59°, and the semi-major axis in units of the stellar 
radius a/Rg = 5.94. As reported in Narita et al. (2009a, 
2009b), these assumptions might lead a certain level of 
systematic errors in results due to uncertainties in the 
fixed parameters, especially in i and a/ Rs. We estimated 
such systematic errors in section 4. We also fixed the 
transit ephemeris = 2454230.9053 [HJD] and the or- 
bital period P = 3.553945 days based on S09. Although 
this ephemeris had an uncertainty of 3 minutes for the ob- 
served transit, the uncertainty was well within our time- 
resolution (exposure time of 12-30 minutes and readout 
time of 1 minute) and thus negligible. The adopted pa- 
rameters above are summarized in table 2. 

Our model had 3 free parameters describing the TrES-4 
system: the radial velocity semiamplitude (K), the sky- 
projected stellar rotational velocity (Vsinls), and the sky- 
projected angle between the stellar spin axis and the plan- 
etary orbital axis (A). We also added two offset velocity 
parameters for respective radial velocity datasets (vi : our 
Subaru dataset, V2. Keck in M07). Note that we fixed 
the eccentricity (e) to zero at first, and the argument of 



periastron (w) was not considered. The assumption of 
zero eccentricity is reasonable since Knutson et al. (2009) 
constrained ecosw < 0.0058 (3(t) based on Spitzcr obser- 
vations of the secondary eclipse. 
We then calculated the statistic 



X 



E 



" ^z.calc 



(2) 



where Wi^obs were the observed radial velocity data and 
I'i.caic were the values calculated based on a Keplerian 
motion and on the RM formula given above, ai were cal- 
culated by the quadrature sum of the internal errors of the 
observed radial velocities and expected stellar jitter level 
of 4.4 m s^^. We note that we used the mean value of the 
case ofAMv < 1, AFc,„ <0.6,5-l^<0.6 in Wright (2005) 
for the jitter, based on the stellar properties reported by 
S09. We determined optimal orbital parameters by min- 
imizing the statistic using the AMOEBA algorithm 
(Press et al. 1992). We estimated la uncertainty of each 
free parameter based on the criterion Ax^ = 1.0. 

4. Results 

We first fitted all Subaru radial velocity samples with 
the published 4 radial velocities presented in M07. The 
upper panels of figure 2 show the radial velocities as a 
function of orbital phase plotted with the best-fit model 
curve, and the lower panels plot the radial velocities as a 
function of HJD. Best-fit parameter values, uncertainties, 
and the reduced for this fit are summarized in the left 
column of table 3. The fit is poor, with = 47.4 and 
30 degrees of freedom {xt = 1-58). Apparently there is 
an inconsistency between the Subaru out-of-transit veloc- 
ities, and the M07 Keck velocities, even after allowing for a 
constant offset between these data sets. Comparison of all 
the observations suggests that there are additional out-of- 
transit radial velocity variations of about 20 m s~^, which 
is the root-mean-squared (rms) residual of the Subaru and 
Keck data. Although the reason for this excess variabil- 
ity is still unclear, some possible reasons are long-term in- 
strumental instabilities, starspots, intrinsic motions of the 
stellar surface ("stellar jitter"), and the presence of other 
bodies in the TrES-4 system. We discuss these possibil- 
ities in section 5. Here, we concentrate on a reasonable 
model of the RM effect in spite of the excess variability, 
to give our best estimate for the sky-projected spin-orbit 
angle of the transiting planet TrES-4b. 

Our approach was to use the Subaru data from the 
transit night (UT 2007 July 13), along with the Keck 
data of M07. We chose not to fit the Subaru data ob- 
tained on the other 3 nights outside of transits. The left 
panel of figure 3 plots the data and the best-fit curve 
(the solid line), and the right panel shows a zoom of the 
RM effect. The results are summarized in the middle col- 
umn of table 3. In this case, the fit is good {xt ~ 0.45) 
and the best-fit model indicates a good spin-orbit align- 
ment, with A = 7.3° ±4.6°. The best-fit model also gives 
Vsin/s = 8.3 ±1.1 km s^^, in agreement with the S09 



4 



N. Narita et al. 



result based on a spectroscopic line analysis. We note 
that some of the early RM analyses using the analytic 
formula for the RM effect (e.g., Ohta et al. 2005; Gimenez 
2006) tended to overestimate the stellar rotational veloc- 
ity ysin/s (e.g., Winn et al. 2005). Recently Hirano et al. 
(2010) and Collier Cameron et al. (2009) studied the rea- 
son for the discrepancy, and Hirano et al. (2010) reported 
an improved method to address this problem. The agree- 
ment between our RM results and the spectroscopic anal- 
ysis of S09 suggests that the RM calibration process ex- 
plained above works well. As for the radial velocity semi- 
amplitude, the fit indicates K = 94.9 ± 7.2 m s"^ (the M07 
result was X = 97.4 ± 7.2 m s^^ for reference). 

In addition, we estimated the sensitivity of our results 
to the choices of the fixed photometric parameters by 
fitting the radial velocities using other choices for those 
parameters: (1) a/Rg = 6.15, i = 82.99° (correspond- 
ing to 1(7 lower limit of the impact parameter in S09); 
(2) a/Rs = 5.73, i = 82.19° (corresponding to Icr upper 
limit of the impact parameter in S09). Consequently, 
we found that respective results for A and V sin Is arc; 
(1) A = 7.2° ±4.5° and F sin/, = 8.6 ± 1.1 km s^^; (2) 
A = 7.5° ±4.8° and Fsin/, =8.0 ±1.1 km s"^ Thus the 
conclusion of the spin-orbit alignment does not change, 
and the uncertainties in the photometric parameters have 
very little effects on our results. 

We noticed that the residuals of the Subaru dataset 
have a small positive gradient with time (see the right 
middle panel of figure 3). This indicates that the Subaru 
data alone prefer a slightly smaller value of K than the 
joint fit of Subaru and Keck data. We therefore fit- 
ted the radial velocities of the Subaru UT 2007 July 13 
dataset only, for reference. The best-fit model is shown 
by the dotted lines in figure 3, and derived parameters 
are summarized in the right column of table 3. The 
data indicate a smaller radial velocity semiamplitude, but 
with a larger uncertainty: K = 64.6 ±27.7 m s~^. This 
result may be supporting evidence that a true K for 
the TrES-4b is actually smaller, although it is not very 
convincing. In this light, it is interesting to point out 
that TrES-4b was previously known as the lowest density 
planet (p = 0.202;° °i^^ g cm'^: S09). Recently WASP- 
17b (p = 0.0921°;°^^ g cm-^: Anderson et al. 2010) and 

Kepler-7b {p = 0.166l:J];[j2o g cm'^: Latham et al. 2010) 
were reported to have lower densities than TrES-4b, and 
thus TrES-4b is currently the third lowest density planet 
discovered so far. The observed out-of-transit radial veloc- 
ity variation and the radial velocity gradient around the 
transit phase may suggest a lower density of TrES-4b than 
previously reported. Thus it is important to obtain radial 
velocities of TrES-4 not only to confirm the additional ra- 
dial velocity variation, but also to measure the density 
of TrES-4b more precisely. We note that this fit gives a 
small reduced chi-squarcd {x1 — 0.35) and A = 5.3° ±4.7°, 
indicating a good spin-orbit alignment as before. 

We adopt a compromise value of the latter 2 models 
as our representative result, namely A = 6.3° ± 4.7°. The 
small difference between the results of the latter 2 models 



shows that there is a small systematic uncertainty due 
to the choice of model. Consequently, we conclude that 
the TrES-4b has a small value of A, based on the model 
of the RM effect. However, since we could not find a 
satisfactory solution that explains the all observed radial 
velocities (including all of the Subaru data and the M07 
data), further radial velocity measurements are desired 
to understand any excess variability and to give greater 
confidence to the orbital solution. 

5. Discussion 

5.1. Possible Causes of the Additional Radial Velocity 
Variation 

Since we could not find an appropriate model for the 
all observed radial velocities at this point (see figure 2), 
we here discuss possibilities of systematic effects as well 
as real sources of excess radial velocity variation. 

- Instrumental Instability of the Subaru HDS 

Since the Subaru radial velocities were gathered on a 
few days in clusters spanning about 1 year, it is of utmost 
importance to know the instrumental stability of the HDS. 
For observations within a single night, Narita et al. (2007) 
studied the radial velocity standard star HD 185144 and 
found that the Subaru HDS is stable within a few m s~^. 
In addition, Winn et al. (2009a) reported that the Subaru 
HDS is stable within a few m s~^ over two weeks based 
on HAT-P-7 observations. Likewise, Johnson et al. (2008) 
did not find systematic offsets for the HAT-P-1 system 
over approximately 1 month, using the same setup (and 
even some of the same nights) as the TrES-4 observations 
presented here. However, specific stability over 1 year 
has not yet been confirmed through monitoring of radial 
velocity standard stars, although studies for such long- 
term stability of the Subaru HDS are in progress. Thus 
we note that the instrumental instability of the Subaru 
HDS is one of the prime possibilities of a cause of the 
additional radial velocity variation at this point. 

- Starspots 

One possibility involves starspots on the photosphere 
of TrES-4. Since the rotational velocity of TrES-4 is rel- 
atively fast (V sin Is = 8.5 km s~^: S09), stellar spots of 
similar size to a planet would cause an apparent radial ve- 
locity shift like the RM effect, on a timescale of the stellar 
rotation period (Prot ~ 11 days, assuming sin/, « 1). For 
example, a dark spot of approximately the same size as the 
planet would lead to a maximum shift of 85 m s~^, while 
smaller spots with less contrast would contribute smaller 
velocities. If this is the case, all the RV data are affected. 
However, one would not expect a hot F8 star to have large 
spots. It is because M07 did not report such possibility of 
stellar spots from the TrES transit survey, S09 reported 
no active Ca HK line emission {logR'-^-^ = —5.11 ±0.15), 
and Knutson et al. (2009) concluded that spot activity 
is unlikely (but not impossible) based on Spitzer obser- 
vations. Thus the spot explanation is doubtful, although 
further long-term photometric monitoring would be useful 
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to constrain this hypothesis still further. 

- Other Sources of Stellar Jitter 

Another possible explanation of the systematic radial 
velocity variation is an intrinsic stellar jitter (see e.g., 
Wright 2005), i.e., motions of the stellar photosphere due 
to pulsations or other flows. Although an empirical rela- 
tion reported by Wright (2005) predicted a typical stellar 
jitter of 4.4 m s~^ for stars like TrES-4, it is conceivable 
that TrES-4 has an unusually unstable photosphere. To 
explain the observed radial velocities, we would need to 
invoke a stellar jitter of about 20 m s~^ for TrES-4 based 
on the rms residuals of the Subaru and Keck datasets. 
The jitter would need to have a time scale longer than a 
few days, in order to explain why the MOT observations 
(conducted on consecutive three nights) do not exhibit 
such a large scatter. In this light the hypothesis that the 
stellar jitter explains all the excess variability seems too 
contrived. 

- Contamination of Companion Star's Lights or Sky 
Backgrounds 

As described by Butler et al. (1996), a slight change in 
the instrumental profile would result in a systematic shift 
of the apparent radial velocity. Thus, any contaminating 
light from the nearby candidate companion star of TrES- 
4 or sky backgrounds (e.g., moonlight or twilight) may 
have affected the radial velocities. We limited the aper- 
ture width of echelle orders in order to avoid contamina- 
tion of lights from the north companion star. Moonlight 
contamination is unlikely since our TrES-4 observations 
were conducted in clear and moonless time, and sky back- 
ground levels were still low although 2 exposures (HJD of 
2454317.74338 and 2454535.14401) were conducted during 
twilight. Note that M07 did not report the existence of 
the companion star, and therefore it is possible that the 
companion star might have been on the slit during the 
M07 observations. Such an effect might have caused some 
systematic shifts in the M07 data. Although we were not 
able to estimate the systematic effect in the M07 data, it 
could be a small effect since the companion is very faint. 

- Eccentricity of TrES-^b 

As M07 reported only 4 radial velocity samples, they 
did not include the eccentricity e and the argument of pe- 
riastron lo in their radial velocity model. Instead they as- 
sumed the orbit to be circular as we have done, and more 
recently Knutson et al. (2009) found ecosw < 0.0058 (3cr) 
based on Spitzer secondary eclipse observations. With the 
Subaru data we now have a sufficient number of radial ve- 
locity samples to allow the eccentricity and argument of 
pericenter to be free parameters. However, allowing e and 
w to vary does not improve the model fit, and the eccen- 
tricity of the best-fit model is nearly zero. This is consis- 
tent with the constraint by Knutson et al. (2009). Thus a 
large eccentricity of TrES-4b could not be the explanation 
for the observed excess RV variability. 



- Additional Planets 

If the preceding explanations for the observed radial ve- 
locity variation could be ruled out, we would consider a 
possibility of presence of additional planets. This is the 
most interesting case, however, at this point our Subaru 
observations were too sparse to find and confirm another 
periodicity in the radial velocities. We can only conclude 
that the radial velocity semiamplitude of about 20 m s~^ 
over a year is possible for hypothetical planets in the 
TrES-4 system. Obviously, further continuous radial ve- 
locity monitoring would be necessary to check on this pos- 
sibility. 

6. Summary 

We observed radial velocities of TrES-4 including a full 
transit of TrES-4b with the HDS of the Subaru 8.2m tele- 
scope. Our radial velocity modeling has revealed that 
TrES-4b has a close alignment between the projected or- 
bital and stellar spin axes, based on the RM effect. On 
the other hand, we could not find a satisfactory Keplerian 
model that agrees with all the data. Although the true 
cause of the excess radial velocity variation is still unclear, 
systematic errors in long-term Subaru observations, stel- 
lar spots, a large stellar jitter, or additional planets might 
play a role in the observed radial velocities. The small 
spin-orbit alignment angle as well as the small eccentric- 
ity of TrES-4b seems to match migration models consid- 
ering disk-planet interactions rather than planet-planet 
scattering models or Kozai migration models. Although 
TrES-4 has the binary candidates to the north and the 
west-southwest, we did not find supporting evidence of 
the Kozai migration. To confirm and discriminate the true 
cause of the radial velocity variation of TrES-4, further ra- 
dial velocity monitoring and photometric monitoring are 
highly desired. 
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the National Astronomical Observatory of Japan. N.N. 
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Table 1. Radial velocities obtained with the Subaru/HDS. 



Time [HJD] 


Value [m s ^] 


Error [m s "'"] 


2454294.7715 


-5.1 


18.6 


2454294.7901 


-12.7 


16.2 


2454294.7991 


-7.3 


16.7 


2454294.8081 


4.5 


15.8 


2454294.8171 


44.1 


16.3 


2454294.8261 


22.5 


17.8 


2454294.8351 


20.7 


15.4 


2454294.8441 


26.8 


16.1 


2454294.8531 


-0.2 


16.6 


2454294.8621 


-10.2 


13.7 


2454294.8711 


-20.9 


14.6 


2454294.8801 


-39.3 


13.8 


2454294.8891 


-52.5 


12.8 


2454294.8981 


-51.6 


14.5 


2454294.9071 


-45.5 


15.2 


2454294.9161 


-77.1 


14.6 


2454294.9252 


-70.0 


12.2 


2454294.9342 


-43.9 


13.0 


2454294.9432 


-41.3 


11.9 


2454294.9522 


-16.5 


14.0 


2454295.0289 


-31.6 


14.2 


2454295.0379 


-36.5 


13.1 


2454295.0559 


-33.5 


13.6 


2454317.7434 


-38.4 


15.8 


2454317.9553 


12.6 


19.6 


2454535.1010 


-10.9 


19.3 


2454535.1225 


-28.2 


17.3 


2454535.1440 


-18.4 


18.0 


2454616.8015 


3.0 


18.8 


2454616.8162 


-24.2 


19.6 


2454616.8309 


-16.4 


20.3 



Table 2. Adopted stellar and planetary parameters. 



Parameter 


Value 


Source 


[Mq] 


1.404 


S09 


Rs [Rq] 


1.846 


S09 


Rpl Rs 


0.09921 


S09 


^n 


82.59 


S09 


a/Rs 


5.94 


S09 


Ul 


0.46 


Claret (2004) 


U2 


0.31 


Claret (2004) 


jitter [m s^-^] 


4.4 


Wright (2005) 


Tc [HJD] 


2454230.9053 


S09 


P [days] 


3.553945 


S09 




Fig. 1. HDS slit viewer images of TrES-4 taken on UT 2007 July 13 (upper left), UT 2007 August 05 (lower left), UT 2008 March 
9 (upper right), and UT 2009 July 12 (lower right). North is up and east is left, and field of view is 20" X 20" for the all images. 
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Fig. 2. Radial velocities (RVs) and the best-fit curve of TrES-4 as a function of orbital phase (upper) and as a function of HJD 
(lower). All Subaru RVs and the M07 RVs are used. Bottom panels: Residuals from the best-fit curve. 
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Fig. 3. Top panels: Radial velocities (RVs) and best-fit curves of TrES-4 as a function of orbital phase. The left panels show the 
entire orbit and the right panels are the zoom of transit phase. The solid line indicates the best-fit curve for the case that Subaru 
RVs taken on UT 2008 July 13 and the M08 RVs are used, and the dotted line is for the case that only Subaru RVs are used. Middle 
panels: Residuals from the solid model curve. Bottom panels: Residuals from the dotted model curve. 



Table 3. Best-fit values and uncertainties of the free parameters. 





Subaru all + Keck 


Subaru transit + Keck 


Subaru transit only 


Parameter 


Value 


Uncertainty 


Value 


Uncertainty 


Value 


Uncertainty 


K [m s"^] 


75.5 


±6.3 


94.9 


±7.2 


64.6 


±27.7 


Vsinlg [km s^^] 


8.3 


±1.1 


8.3 


±1.1 


8.7 


±1.2 


A[°] 


0.0 


±4.2 


7.3 


±4.6 


5.3 


±4.7 


vi (Subaru) [m s ^] 


-23.2 


±3.8 


-13.0 


±4.6 


-15.5 


±5.1 


rms (Subaru) [m s^^] 


19.9 




9.2 




8.4 




V2 (Keck) [m s^^] 


18.7 


±5.9 


22.9 


±6.0 






rms (Keck) [m s^"'^] 


18.4 




7.9 










47.36/30 




9.88/22 




6.64/19 





